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ABSTRACT 


Increasing  awareness  of  the  applicability  of  systems  analysis 
to  urban  problems  results  in  a program  to  develop  concepts  which  are 
relevant  to  locational  problems  in  various  urban  areas.  The  goal  of 
this  project  is  the  development  of  a model  to  provide  the  urban  decision 
maker  with  quantitative  information  regarding  alternative  choices  in 
determining  public  service  facility  locations.  In  the  case  of  fire 
station  location,  this  model  evaluates  possible  location  sites  by  testing 
sensitivity  to  demographic  and  structural  developments  in  the  area 
served  by  the  facilities. 
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A Computer  Model  for  Evaluating  Fire  Station  Location 

I . Summary 

In  January  1968,  the  International  City  Managers’  Association 
(ICMA)  sponsored  a symposium  to  discuss  the  application  of  systems 
analysis  to  the  resolution  of  urban  problems.  The  symposium  generated 
interest  in  a program  to  evaluate  the  applicability  of  systems  analysis 
to  specific  problems.  Program  sponsorship  was  under  the  auspices  of  the 
ICMA,  the  American  Society  of  Planning  Officials,  the  Fels  Institute  of 
Local  and  State  Government  of  the  University  of  Pennsylvania,  and  the 
Technical  Analysis  Division  of  the  National  Bureau  of  Standards. 

This  report  is  addressed  to  one  project  done  within  the  constraints 
of  the  program.  This  project  provided  a means  for  analyzing  the  location 
of  public  service  facilities.  East  Lansing,  Michigan,  was  selected  as 
the  participant  city  for  study  and  the  Technical  Analysis  Division  was 
designated  to  serve  as  "technical  coach"  for  the  work  to  be  done.  The 
goal  of  the  project  was  to  generate  information  that  could  be  used  by 
the  urban  decision  maker  in  the  process  of  locating  fire  stations  according 
to  present  and  future  protection  requirements  of  the  city. 

In  the  case  study  work  done  for  East  Lansing,  the  specific  objectives 
were  defined  as: 

‘to  evaluate  fire  protection  needs  of  the  city  in  the  static  situation 
of  the  present  environment  and  to  determine  the  optimal  number  and 
location  of  fire  stations  requisite  to  meet  these  needs; 


•to  determine  the  protection  needs  of  the  city  as  of  a future  date 
(1980)  and  to  ensure  that  fire  stations  which  are  adequate  today 
will  continue  to  meet  the  protection  requirements  of  the  changing 
environment . 

The  city  staff  of  East  Lansing  and  TAD  analysts  cooperated  in  the 
development  of  a simulation  model  to  accomplish  these  objectives.  In 
the  process,  it  was  demonstrated  that  systems  analysis  can  be  applied  to 
a city  problem  at  a reasonable  cost  and  that  the  analysis  can  be  performed 
by  the  city  staff.  It  is  hoped  that  city  staffs  in  other  regions  of  the 
country  who  face  similar  public  service  facility  location  problems  will 
be  able  to  utilize  this  model  with  little  modification. 

As  a result  of  the  work  done  in  East  Lansing,  the  coverage  area  of 
a fire  station,  located  on  the  Michigan  State  University  campus,  is  being 
expanded;  a second  station  is  being  relocated  north  of  the  central  business 
district;  and  a new  station  will  be  constructed  in  the  Southeast  Section 
of  the  city  (see  figure  3) . It  was  determined  that  these  three  fire 
stations  will  afford  adequate  protection  for  the  changing  city,  now  and 
during  the  coming  decade. 

The  procedure  in  East  Lansing  consisted  of  analyzing  trial  locations 
known  to  be  reasonable,  based  on  land  availability,  land  cost,  and  citizen 
consent.  A systematic  analysis  of  the  results  met  the  prime  objectives 
of  the  study. 
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II.  THE  MODEL 


The  model  is  designed  to  evaluate  the  effectiveness  of  an  existing 
or  trial  station  location.  Mathematically,  the  model  is  structured  around 
a shortest-path  algorithm,  which  calculates  response  time,  and  a "weighting 
function,"  which  quantifies  or  measures  the  fire  hazard  of  each  structure. 

The  various  members  of  the  Project  Team  agreed  that  response  time  and 
fire  hazard  are  the  two  most  relevant  factors  in  determining  station 
location . 

The  present  criteria  for  locating  fire  stations  are  determined  by 
the  American  Insurance  Association.  The  distance  criterion  establishes 
maximum  limits  in  terms  of  travel  distance  to  either  high  value  or 
residential  areas.— ^ This  distance  criterion  is  most  applicable  to 
those  cities  where  high  value  construction  is  in  the  central  business 
district  surrounded  by  a symmetric  development  of  residential  homes. 

This  is  an  unrealistic  view  of  the  structural  and  demographic 
characteristics  of  cities  today.  Also,  present  transportation  systems 
do  not  lend  themselves  to  a simple  correlation  between  distance  and  time. 
That  is  to  say,  a one -mile  drive  on  an  expressway  obviously  requires  less 
travel  time  than  a one-mile  drive  in  a central  business  district.  There 
is  need  for  the  criteria  to  be  sensitive  to  today's  environmental  conditions. 


— ^"The  Grading  of  Municipal  Fire  Protection  Facilities,"  The  National 
League  of  Cities. 
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A.  Shortest-Path  Algorithm 


To  utilize  the  shortest-path  algorithm,  first  it  is  necessary  to 
develop  an  adequate  representation  of  the  street  network  (that  configuration 
of  streets  comprising  the  given  urban  area) . The  network  is  represented 
by  a uniquely  numbered  series  of  nodes  and  links.  A node  indicates  each 
street  intersection  worthy  of  representation,  i.e.,  it  is  not  necessary 
to  represent  small  street  intersections.  A link  represents  the  distance 
and  travel  time  between  each  pair  of  nodes. 

In  many  parts  of  the  country,  the  city,  county,  or  state  highway 
department  has  coded  the  street  network  for  present  and  future  years. 

This  representation  is  usually  prepared  according  to  a Bureau  of  Public 
Roads  format.  The  data  listed  include  distance  between  intersections 
(links  between  nodes),  off-peak,  peak,  and  posted  speed  to  traverse  each 
link,  width  of  roadway,  etc. 

In  the  East  Lansing  study,  the  Tri -County  Highway  Department  had 
coded  the  East  Lansing  street  network  for  1965  and  1975. 

Having  adequate  network  representation,  the  algorithm  was  used  to 
calculate  the  minimum  time  path  from  one  node , representing  a trial  fire 
station  location,  to  all  other  nodes.  The  East  Lansing  Fire  Department 
indicated  that  the  travel  time  for  their  engines,  between  any  pair  of 
nodes,  was  well  approximated  by  dividing  the  distance  by  the  posted  speed, 
and  neglecting  any  change  of  speed  at  an  intersection.  However,  there 


-4- 


are  two  continually  congested  intersections  in  East  Lansing  and  a 
half -minute  of  travel  time  was  added  to  all  trips  traversing  either 
of  these  intersections. 

B.  The  "Weighting  Function" 

Appraisal  of  the  shortest-path  times  affords  one  means  for  evaluating 
potential  fire  station  locations.  However,  some  properties  are  more 
important  than  others  in  terms  of  potential  damage.  This  variation  in 
fire  hazard  is  measured  by  a "weighting  function."  The  "weighting 
function"  may  be  defined  as  a heuristic  function  which  quantifies  a 
judgment  of  the  East  Lansing  Project  Team.  This  judgment  concerns  the 
relative  significance  of  each  of  the  structural  and  human  factors  as 
related  to  the  hazard  of  a potential  fire  on  a given  property. 

The  rationale  for  such  a "weighting  function"  is  clearly  seen 
when  considering  the  protection  requirements  of  an  empty  lot  vs.  a lot 
containing  a school.  Thus,  the  East  Lansing  Project  Team  decided  that 
the  "weighting  function"  should  be  a combination  of  the  following 
factors : 

. the  land  use  type; 

. the  potential  danger  to  human  life;  and 

. the  potential  loss  of  the  structure. 

These  factors  contain  a number  of  variables  which  had  to  be  identified, 
their  interrelationships  defined,  and  their  value  quantified.  The 
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Project  Team  identified  the  critical  variables  as: 

1)  number  of  persons  endangered; 

2)  number  of  stories  in  each  building; 

3)  construction  type; 

4)  age  of  structure;  and 

5)  square  area  of  each  floor. 

The  "weighting  function”  was  devised  by  systematically  correlating 
structural  and  demographic  characteristics  of  all  land  use  types  with 
the  statistical  results  of  an  analysis  of  the  city’s  fire  history  data. 

The  explicit  physical  characteristics  of  each  structure,  as  well  as 
implicit  characteristics  according  to  land  use  classification,  were 
quantified  as  they  related  to  the  need  for  fire  protection.  The 
final  product,  (Eq.  1)  provides  a definite  correlation  between  the 
measurable,  structural  and  demographic  characteristics  and  frequency 
of  fire.  The  frequency  of  fire  is  the  percentage  of  fire  by  occupancy 
type  in  proportion  to  the  total  number  of  structures  of  that  land  use 
type.  The  close  correlation  between  the  frequency  of  fire  and  the 
total  points  assessed  for  the  various  occupancy  types  experimentally 
suggests  the  validity  of  the  penalty  point  system  as  a quantitative 
description  of  the  fire  hazard  (see  Chart  1) . The  formulation  of  the 
"weighting  function"  is: 

Structure  = population  (1)  + height  in  stories  (2)  + 

(Eq.  1) 

construction  type  (3)  + age  (4)  + area  in  square  feet  (5) 

Weighting  function  = structure  X occupancy  type  (6).  (Eq.  2) 
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Chart  1 - The  Weighting  Function 
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2/  Churches  are  being  used  more  as  schools,  and  therefore,  given  the  same  land  use  value. 


(1)  Population  is  the  number  of  persons  who  actually  occupy  the 
dwelling.  Apartments  and  single-family  dwelling  units  have  been  given 
an  average  population  per  household  of  2.5  and  3,  respectively. 

(2)  Height  in  stories  has  been  penalized  as  follows:  One  and  two 

stories,  0 points;  three  and  four  stories,  2 points;  five  and  six 
stories,  4 points;  seven  and  over,  10  points. 

(3)  Construction  type  of  each  structure  is  placed  on  a weighted 
scale  with  a range  of  four  penalty  points.  The  critical  end  of  the 
scale  is  4,  and  all  single-family  dwelling  units  have  been  assigned  a 
M4"  rating. 

(4)  Age  of  structure  has  been  penalized  as  follows:  0 to  10  years, 

1 point;  11  to  20  years,  2 points;  21  to  30  years,  3 points;  31  to  40 

years,  4 points  . . . 

(5)  Area  in  square  feet  has  been  multiplied  by  1/1000  in  order  to 
bring  the  extremely  large  numbers  into  perspective  with  the  other  weighted 

categories . 

(6)  Occupancy  types,  the  multiplicative  factor  in  Eq.  2,  have  a 
value  of  one  (1.0),  except  for  schools  and  churches,  which  have  been 
assigned  the  value  of  one  third  (.333). 

The  link  between  the  physical  description  of  the  city  and  its  socio- 
economic description  is  achieved  by  multiplying  the  value  of  the  "weighting 
function"  for  a structure  by  the  length  of  time  it  takes  an  emergency 
vehicle  to  reach  that  structure.  (See  Appendix  B for  mathematical  formu- 
lation.) 
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East  Lansing  provided  the  Project  Team  with  data  on  all  actual 
and  potential  fires  within  the  city  for  the  past  10  years.  These 
records  contain  the  node  number  location  of  each  fire,  the  date,  day 
of  the  week,  hour,  responding  station,  response  distance,  fire  duration, 
number  of  persons  endangered,  cause,  structure  type,  occupancy  type, 
height  of  building,  year  building  constructed,  fire  rating  of  structure, 
amount  of  damage  to  structure,  amount  of  damage  to  contents,  and 
address  of  structure  (see  Appendix  A) . 

C.  Districting 

The  designation  of  districts,  or  station  coverage  areas,  is 
presently  a manual  operation.  To  evaluate  a trial  location,  it  is 
necessary  to  input  the  coverage  area  for  the  location.  Recognizing 
that  the  area  of  a district  is  reflected  in  the  model  evaluation,  it 
is  necessary  to  reassess  the  district  after  each  effectiveness  rating. 

In  the  East  Lansing  application,  isochronal  boundaries  of  2 1/2  to 
3 minutes  were  used  to  outline  the  first  district  for  each  trial  location. 
Thus,  the  Project  Team  studied  various  location  plans  through  the 
manipulation  of  districts. 

D.  Assumptions 

A basic  assumption  in  this  model  is  that  all  factors  influencing 
the  location  of  a fire  station  can  be  expressed  in  terms  of  a common 
utility  scale.  For  example,  in  simulating  the  street  network,  probabilistic 
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measurement  of  variables  such  as  draw-bridges  and  railroad  crossings 
must  be  expressed  through  the  impedance  value  associated  with  the 
appropriate  link.  However,  the  way  in  which  these  points  are  introduced 
into  the  evaluation  remains  flexible. 

Time  is  used  as  the  metric  for  analyzing  the  highway  network.  The 
time  associated  with  each  link  is  calculated  by  dividing  the  distance 
in  miles  by  the  posted  speed  in  miles  per  hour.  The  travel  time  to  each 
destination  is  the  sum  of  the  travel  times  for  those  links  comprising 
the  shortest  time  path.  It  is  assumed  that  any  speed  variation  within 
a block  or  at  an  intersection  approximates  the  posted  speed.  This 
assumption  was  validated  in  East  Lansing  by  matching  the  computer- 
calculated  time  with  the  actual  response  times  to  fires. 

A third  assumption  concerns  the  aggregation  of  the  "weighting 
functions"  for  each  land  use  to  the  node  in  the  lower  right-hand  comer 
of  each  neighborhood.  Such  an  aggregation  is  necessary  to  solve  the 
problem  with  the  existing  computer  hardware.  It  was  the  opinion  of  the 
East  Lansing  Project  Team  that  such  an  aggregation  did  not  significantly 
affect  the  sensitivity  of  the  model. 

E . Evaluation 

The  effectiveness  measure  for  each  possible  site  is  achieved  by 
taking  the  product  of  response  time  and  "weighting  function"  and  sunming 
it  for  all  points  in  the  district.  By  comparing  the  sum  total  for 
each  trial  site,  a ranking  of  locations  is  obtained. 
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The  final  selection  of  a location  site  is  based  on  several  criteria 


1)  a measure  of  effectiveness  which  is  not  significantly  affected 
by  proposed  development  plans; 

2)  a districting  plan  which  minimizes  excessive  overlapping  of 
coverage  areas  as  additional  stations  are  required;  and 

3)  a satisfactory  trade-off  between  an  acceptable  level  of 
effectiveness  and  the  financial  and  political  constraints. 

By  satisfying  these  three  conditions,  the  city  provides  both  effective 
and  efficient  fire  protection. 
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III.  THE  PROGRAM 


The  basic  input  for  the  program  is  a highway  network  (see  Appendix  C) 
and  a structural  data  index  (see  Appendix  D) . The  network  consists  of 
nodes  and  links.  Each  node  corresponds  to  an  arbitrarily  numbered  inter- 
section on  a map  of  the  city.  In  the  central  business  district,  every 
intersection  is  coded;  while  in  the  less  densely  populated  residential 
area,  every  second  or  third  intersection  is  coded.  Each  link  represents 
a street  segment  between  two  intersections  and  is  identified  by  the  two 

nodes.  The  link  is  characterized  by  travel  direction  (one-way  or  two-way) 

2/ 

posted  speed,  and  distance.  Travel  time  is  calculated  accordingly  to 

distance  and  posted  speed. 

The  form  of  the  data  is  expanded  by  a pre-processing  routine,  so 
that  each  link  is  identified  by  a unique  pair  of  ordered  nodes.  This 
ordering  makes  it  possible  for  travel  time  in  one  direction  to  differ  from 
travel  time  in  the  opposite  direction.  Since  fire  equipment  will  not 
normally  travel  against  traffic,  it  is  important  to  identify  one-way  streets. 
This  pre-processing  routine  also  sequentially  renumbers  the  expanded  list 
of  links,  which  is  necessary  for  the  computer  analysis  of  the  network. 

These  node  numbers  are  referred  to  hereafter  as  the  renumbered  node  numbers . 

77 

— Other  characteristics,  such  as  width  of  roadway,  speed  at  peak  and 
off-peak  hours,  etc.,  are  available  from  an  "Inventory  Summary  for 
Selected  Streets  and  Highways  in  the  Tri -County  Region,"  prepared  by 
the  Lansing,  Michigan  Tri-County  Planning  Organization.  These 
characteristics  were  not  used  in  this  analysis. 
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The  building  index  is  a listing  of  all  land  use  classifications,  i.e., 
buildings  and  vacant  lots  in  East  Lansing.  Each  structure  is  characterized 
by  height,  construction  type,  structural  age,  and  square  footage  per 
floor.  Each  building  in  a neighborhood  defined  by  the  intersection  of 
three  or  more  links  is  associated  with  the  node  in  the  lower  right-hand 
comer  of  that  neighborhood.  These  two  inputs  are  usually  read  from 
magnetic  tape. 

The  data  tape  is  created  by  the  two  programs,  DATAT  and  STRUCT  (see 
Appendix  E) . DATAT  creates  the  network  tape  in  a form  fitted  for  the  input 
requirement  of  program  RANK.  STRUCT  creates  the  structural  index  file 
following  the  network  information.  Consequently,  only  one  tape  is  used 
to  input  information.  The  logical  unit  number  of  the  tape  drive  is  specified 
by  the  variable  TAP.  The  logical  unit  numbers  of  the  card  reader  (for  input) 
and  printer  (for  output)  are  designated  by  the  parameters  R and  P,  respectively. 

The  card  input  which  identifies  trial  sites  and  districts  is  divided 
into  three  sections.  The  first  section  identifies  which  program  options 
are  selected,  sets  the  searching  limits  of  the  shortest -path  algorithm,  and 
initializes  other  boundary  constraints.  The  second  section  specifies  those 
changes  in  structural  development  being  considered;  and  the  third  section 
identifies  the  coverage  area  for  each  set  of  trial  locations. 

A.  Section  1 - Parameters  and  Options 
This  section  contains  eleven  (11)  cards.  The  first  card  identifies 
four  options: 

NRUNS  - the  number  of  passes  through  the  model; 

0PT1  - optional  output  display  of  distance  and  time  from  origin  to 
all  destinations; 
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0PT2  - optional  output  listing  of  all  nodes  and  the  summed  value  of 
the  "weighting  function"  associated  with  each  node;  and 
FUTURE  - option  to  consider  future  development. 

Based  on  a particular  highway  network  and  structural  index  a single  run 
of  the  program  evaluates  a number  of  origins  associated  with  a districting 
plan.  The  variable,  NRUNS,  is  both  the  number  of  passes  through  the  model 
and  the  number  of  districts.  For  each  node  chosen  as  an  origin,  two  arrays 
are  calculated  by  the  shortest-path  algorithm.  In  the  first  array,  each 
element  is  the  shortest-path  travel  time  from  an  origin  to  all  other  destin- 
ations. Each  element  in  the  second  array  gives  the  distance  from  an  origin 
to  all  destinations  when  traveling  along  the  least -time  path.  When  0PT1  is 
triggered,  the  program  will  print  the  subsequent  arrays  for  all  nodes 
selected  as  trial  origins  (see  Section  3).  0PT2  is  both  an  optional  print- 
out of  two  lists  of  nodes  and  the  summed  value  of  the  "weighting  function" 
associated  with  each  node.  One  list  is  the  renumbered  node  numbers  used 
during  the  execution  of  the  program;  the  second  list  contains  the  original 
node  numbers  which  appear  on  the  work  map  of  the  city.  The  fourth  option, 
FUTURE,  signals  a decision  to  consider  future  structural  developments.  When 
this  option  is  triggered,  the  program  will  read  data  concerning  proposed 
structural  development  and  change. 

The  second  card  in  this  section  assigns  values  for  three  variables: 

MAX  - the  limit  on  the  number  of  links  the  algorithm  will  search  to 
determine  a shortest  time  path. 

NNODE  - the  number  of  nodes  in  the  network. 

NLINK  - the  number  of  unique  links  in  the  network. 
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In  the  arithmetic  "weighting  function,"  land  use  is  a multiplicative 
factor.  The  alphanumeric  land  use  classifications  are  contained  in  the 
list  ICODE.  Presently,  the  program  reads  64  land  use  classifications. 

The  numeric  values  assigned  with  each  land  use  type  are  in  the  list,  ASSIGN, 
and  are  referenced  by  the  same  subscript  as  the  corresponding  variable  name. 

The  final  card  in  this  section  defines  the  parameter  values  for  the 
"weighting  function."  These  values  are  contained  in  List  A. 

B.  Section  2 - Future  Development 

This  section  modifies  the  information  from  the  structural  index  tape 
to  incorporate  future  structural  development  and  change.  When  the  option, 
FUTURE,  is  used,  the  analyst  has  the  ability  to  exclude  structures  contained 
on  the  index  list,  and  to  include  additional  structures  that  will  exist  in 
the  future.  The  program  will  delete  any  number  of  structures  of  a particular 
land  use  type  from  the  neighborhood.  At  present,  the  program  can  delete 
four  different  types  of  structures  from  ten  different  neighborhoods.  The 
first  card  in  this  section  gives  the  general  list  of  land  use  classifications 
to  be  removed  as  future  development  occurs.  The  variable  name  of  this 
list  is  NSTR.  Each  of  the  10  cards  following  this  one  identifies  the 
original  node  number  from  which  structures  are  to  be  removed,  and  the 
number  of  structures  to  be  removed  of  each  land  use  type.  The  variable, 
DETROY,  identifies  the  list  containing  the  nodes  from  which  structures  are 
to  be  removed.  LIMIT  is  the  list  containing  the  number  of  structures  to 
be  removed. 
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Following  the  11  cards  described  previously,  the  deck  of  cards 
containing  the  proposed  structures  completes  the  data  (see  Appendix  D) . 

Any  number  of  cards  describing  future  development  will  be  read  at  this 
time.  At  the  end  of  the  deck,  there  must  be  End-of-File  card.  If  future 
development  is  not  to  be  considered  in  the  analysis,  this  section  of  the 
program  is  omitted  entirely. 

C.  Section  3 - Service  Districts 
The  third  section  is  concerned  with  the  identification  of  service 
areas.  The  information  in  this  section  can  be  repeated  a number  of  times 
(NRUNS  times).  The  first  card  is  the  title  card,  identifying  the  area 
under  consideration.  The  next  card  concerns  identification  of  the 
trial  locations.  The  selection  of  nodes  to  be  evaluated  as  test  sites  is 
contained  in  the  list  ICK.  An  option  exists  to  printout  the  series 
of  links  comprising  the  shortest  time  path.  The  user  selects  any  number 
of  nodes  from  the  subset  of  nodes  chosen  as  test  origins.  This  set  of  nodes 
is  contained  in  the  list  ICHO.  The  program  will  enumerate  the  paths  from 
the  selected  node  to  all  other  nodes  in  the  city  designated  in  the  list  IDEST. 
This  information  is  useful  for  validating  the  highway  network. 

In  this  section,  the  first  option,  JSS,  allows  the  user  to  select  a 
specific  area  of  coverage.  Variations  in  existing  equipment  and  manpower 
are  reflected  in  the  analysis  by  selecting  districts  which  are  proportional 
to  available  equipment  and  manpower.  If  JSS  is  activated,  it  is  followed 
by  cards  designating  all  nodes  within  the  service  area.  The  list  containing 
these  selected  nodes  is  called  JSID. 
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When  evaluating  possible  locations  for  a fire  station  having  prime 
responsibility  for  a specific  district,  it  may  be  desirable  to  include  some 
fractional  proportion  of  responsibility  for  protecting  an  area  not  encompassed 
in  the  home  district.  This  option  to  consider  secondary  coverage  is  triggered 
by  the  variable  NPOL.  If  this  option  is  desired,  the  percentage  of  secondary 
coverage  is  specified  by  the  variable  POLICY.  Following  triggering  of  this 
variable,  the  nodes  in  the  area  under  secondary  coverage  protection  are 
enumerated  in  the  list  IPOL.  If  NPOL  is  not  triggered,  the  run  terminates. 
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IV.  USER’S  MANUAL 


The  model  is  basically  a single  program  developed  around  the  shortest- 
path  algorithm.  Although  the  program  contains  two  subroutines  (CARTIM  and 
HOUSE),  they  do  not  perform  a single  unique  operation;  and  consequently, 
the  discussion  of  the  model  will  refer  to  the  model  as  one  program. 


The  program  has  two  types  of  inputs:  tape  and  card.  The  tape 

contains  the  coded  highway  network  and  the  existing  structural  index  (see 
Appendix  E) . The  card  input  initializes  the  parameters  and  identifies 
various  model  options  and  printout  options. 

The  activation  of  an  option  or  the  selection  of  a node  is  accomplished 
by  placing  a one  in  the  proper  column.  The  consequences  of  not  triggering 
an  option  is  discussed  at  the  same  time  as  the  capabilities  of  the  option. 
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In  the  case  of  node  selection  for  the  lists  JSID  and  IPOL,  a node  is 
designated  by  punching  a one  in  the  column  corresponding  to  the  renumbered 
node  identification.  Thus,  the  program  sequentially  reads  in  columns 
on  (NNODE/80)  cards.  Otherwise,  nodes  are  designated  by  the  number  appearing 
on  the  work  map  of  the  city. 

The  card  input  is  divided  into  three  sections  for  discussion  purposes 
only.  The  "parameters  and  options"  section  identifies  the  program  options 
and  sets  parameter  values.  This  secricn  contains  11  cards. 


No.  of  Cards 


1 


1 


Column  (or  Format)  Description 

1-3  NRUNS  - the  number  of  consecutive 

iterations  of  the  program 
changing  areas  of  coverage 
and/or  policy  decisions. 

6*  I0PT1  - the  option  to  print  time  and 

distance  arrays  associated 
with  each  origin. 

9*  I0PT2  - the  option  to  print  a dic- 

tionary of  all  nodes  along 
with  the  weighting  function 
value  for  each  node. 


12* 


1-5 


6-10 


FUTURE  - the  option  to  consider  future 
structural  development. 

MAX  - the  maximum  number  of  links 
investigated  to  build  a 
path. 

NNODE  - the  number  of  nodes  in  the 
network . 


*Place  a one  in  proper  column  to  trigger  option. 
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No.  of  Cards 


Column  (or  Format) 


Description 


11-15 

NLINK  - 

the  number  of  links  in 
the  network. 

16-17 

TAP  - 

logical  unit  number  of 
data  tape. 

18-19 

R - 

logical  unit  number  of 
the  card  reader. 

20-21 

P - 

logical  unit  number  of 
the  printer. 

4 

16A4 

I CODE  - 

the  list  of  land  use 
categories . 

4 

16F4.2 

ASSIGN  - 

the  numerical  value 
assigned  to  each  land 
use. 

1 

9F6.3 

A - 

the  list  of  coefficients 
for  the  weighting  function. 

The  number  of  cards  in  the  following  two  sections  is  a function  of 
the  choice  of  options  activated  in  section  I. 

Section  II  specifies  the  changes  in  structural  development  that  are 
to  be  considered.  If  there  are  no  changes,  this  section  is  omitted.  When 
additional  structural  development  is  to  be  considered  (triggered  by 
assigning  FUTURE  = 1),  then  include  the  following  cards: 

No.  of  Cards  Column  (or  Format)  Description 

1 3A4  NSTR  - a list  of  the  land  use 

types  to  be  removed  when 
future  development  takes 
place. 

10  16,  IX  DETROY  - the  node  from  which 

structures  are  to  be 
removed  (original  node 
number) . 
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No.  of  Cards  Column  for  Format)  Description 

313  LIMIT  - the  number  of  structures 

of  each  land  use  type 
in  NSTR  to  be  removed 
(start  column  8) . 

Place  future  development  deck  of  cards  here,  followed  by  an  End- of - 
File  card. 


The  "service  districts"  section  identifies  the  set  of  trial  origins 
and  their  respective  coverage  areas  for  each  pass.  The  following  cards 
are  changed  according  to  the  number  of  passes  (NRUNS)  during  a single 
run  of  the  program. 


Cards 

Column  (or  Format) 

Description 

1 

12A6 

TITLE 

- identification  of  the 
coverage  area. 

1 

10(15) 

ICK 

- nodes  chosen  to  be  origins 
(original  node  number) . 

1 

10(4X,  11)* 

ICHO 

- nodes  chosen  from  which  to 
enumerate  paths  to  all 
destinations . 

1 

10(15) 

IDEST 

- nodes  chosen  as  path 

destinations  (original  node 
number) . 

1 

3* 

JSS 

- option  to  consider  a 
specific  area  rather  than 
the  entire  city.  If  zero, 
skip  to  NPOL. 

6 

8011 

JSID 

- nodes  within  the  specified 
area. 

1 

3 

NPOL 

- option  to  consider  an  area 
of  secondary  coverage.  If 
zero,  this  is  the  end. 

1 

F4.2 

POLICY 

- the  percentage  of  secondary 
coverage . 

*Place  a one  in  proper  column  to  trigger  option. 
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No . of  Cards 


Column  (or  Format) 


Description 


6 8011*  IPOL  - nodes  considered  under 

policy  percentage  of 
secondary  coverage. 

For  another  pass,  the  next  card  following  this  section  will  be  another 
title  card.  Of  course,  the  subsequent  cards  would  then  follow. 

The  computer  program  has  been  implemented  in  East  Lansing,  Michigan. 
The  program  is  written  in  FORTRAN  IV  for  the  CDC  3600  at  MSU  and  FORTRAN  V 
for  the  UNIVAC  1108  at  NBS.  Run  on  the  CDC  3600  to  evaluate  a network  of 
432  nodes,  1400  links,  and  1500  structures,  the  program  evaluated  10  trial 
sites  in  each  of  three  districts  in  approximately  3 1/2  minutes. 


*Place  a one  in  proper  column  to  trigger  option. 
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V.  LEXICON 


I . Programs 

DATAT  - computer  program  to  transfer  street  network  to  the  data  tape. 
RANK  - computer  program  which  evaluates  the  effectiveness  of  trial 
fire  locations. 

STRUCT  - computer  program  to  transfer  structural  and  demographic 
information  to  the  data  tape. 

II.  Single  Parameters 

FUTURE  - the  option  to  consider  future  structural  development 
JSS  - the  option  to  consider  a specific  service  area. 

MAX  - the  maximum  number  of  links  investigated  by  the  shortest-path 
algorithm  to  build  a least- time  path. 

NLINK  - the  number  of  unique  links  in  the  street  network. 

NNODE  - the  number  of  nodes  in  the  network. 

NPOL  - the  option  to  consider  an  area  of  secondary  coverage. 

NRUNS  - the  number  of  consecutive  iterations  of  the  program. 

0PT1  - the  option  to  print  time  and  distance  arrays. 

0PT2  - the  option  to  print  a dictionary  of  all  nodes  along  with 
the  weighting  function  value  for  each  node. 

P - the  logical  unit  number  of  the  printer. 

R - the  logical  unit  number  of  the  card  reader. 

TAP  - the  logical  unit  number  of  the  data  tape. 
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III.  Lists 


A 

ASSIGN 

DETROY 

ICHO 

ICK 
I CODE 
IDEST 
I POL 
JSID 
LIMIT 
NSTR 

POLICY 

TITLE 


the  coefficients  for  thw  weighting  function. 

the  numerical  value  assigned  to  each  land  use. 

the  nodes  from  which  structures  are  to  be  removed. 

the  nodes  chosen  from  which  to  enumerate  paths  from  trial 

origins  to  selected  destinations. 

the  nodes  chosen  to  be  origins. 

the  land  use  categories . 

the  nodes  selected  as  destinations. 

the  nodes  under  secondary  coverage. 

the  nodes  within  the  specified  service  area. 

the  number  of  structures  of  each  land  use  type  to  be  removed. 

the  land  use  types  to  be  removed  when  future  development 

takes  place. 

the  percentage  of  secondary  coverage, 
the  identification  of  the  coverage  area. 
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City  of 
EAST  LANSING 


21/2  minute  Ischronal 
Lines 

Present  Fire  Locations 


City  of 
East  Lansing 

* 


Proposed  Boundary  Lines 

Proposed  Fire  Station 
Locations 
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Figure  3 


Appendix  A - Fire  Record  Code  Format 


Column  Information 


2 

- 5 

Node  Number 

8 

- 12 

Date  (8-month;  9 § 10 -day  of  month;  11 
12 -year) 

13 

Day  of  Week 

14 

Hour  (nearest) 

15 

a.m.  or  p.m. 

17 

Station  first  responding 

19 

- 20 

Distance  in  lOths  of  a mile 

23 

- 25 

Duration  of  fire  in  minutes 

27 

- 30 

Number  of  persons  endangered 

32 

Cause 

34 

- 35 

Structure  occupancy  type 

39 

- 40 

Height  of  building  (number  of  stories) 

42 

- 43 

Year  building  constructed 

45 

Fire  rating  of  structure 

49 

- 55 

Amount  of  damage  to  structure 

59 

- 65 

Amount  of  damage  to  contents 

73 

- 80 

Address  (73-77  structure  number;  78-80 

street  code) 


A-l 


Appendix  B - Mathematical  Formulas 


1)  The  weighting  function  ^ for  each  building  k associated  with 


node  j can  be  expressed  as: 

5 


) /"I  • n Z 0 1 • X • s i 

(k,j)  k i=1  i i (k,j) 


where 


2,^  = a factor  depending  on  land  use 

th 

a.  = coefficient  of  l 1 category 

X- = i^  structural  characteristic  of  building  (k/,j) 
1 J 

2)  The  aggregated  value  associated  with  node  j is: 


m 

W . = E .. 

j k=i  (k,n 


where 

n = number  of  buildings  identified  with  node  j 
3)  The  objective  function  F for  origin  3 

p 


f*  ■ L "j 


where 

t . = shortest  path  time  from  node  8 to  node  i. 
8 j 

N = number  of  nodes  in  the  coverage  district. 
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Appendix  C - Highway  Network  Format 


Column 

Information 

1 - 5 

NODE  A 

8 - 12 

NODE  B 

15  - 18 

DISTANCE 

21 

DIRECTION 

24  - 25 

POSTED  SPEED 

1.  NODE  A - the  first  node  used  to  identify  a link  []  D D 0 D 

2.  NODE  B - the  second  node  required  to  identify  a link  jj  P T []  jj 

3.  DISTANCE  - link  distance  in  hundredths  of  a mile  J_|  • Q D 

4.  DIRECTION  - j 0 | - if  one  way  in  the  opposite  direction 

1 1 1 - if  one  way  street 

. 1/ 

] 2 | - if  link  is  a two  way  street  — 


— ^ 3 was  used  instead  of  2 for  two  way  streets  in  the  1975  forecasted 
network. 
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Appendix  D - Structural  Data  Format 


Column 


1 ■ 7 Structural  value  in  thousands  Q Q D D D • [] 


9 - 13 

Land  value  in  thousands  2/  QDD  • D 

18  - 21 

Land  use  □ Q []  □ 

25  - 28 

Structural  age  Q fl  []  Q 

31  - 35 

Square  footage  in  thousands  3/  []  Jj  [] 

39  - 42 

Height  of  building  []["]•  Q 

47  - 49 

Structure  rating  [][][] 

51  - 55 

Number  of  people  Q[]QriDQ 

76  - 80 

Node  number  Q D D D HI 

Example: 

1.  Structural  Value 

Structural  value  will  be  recorded  in  hundreds  and  will,  therefore, 
require  five  columns. 

Examples:  $1,337,800.00 

$ 3,200.00 


1 

3 

3 

7 

8 

0 

0 

0 

3 

2 

2.  Land  Value 


Land  values  will  be  designated  in  hundreds  and  will  require  four  columns 
Example:  $1,400.00  1 |4  | 


2/  Land  values  are  assessed  at  1/3  of  real  value 

3/  The  number  of  square  feet  given  in  each  case  is  only  for  the  ground 
floor 
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3.  Land  Use/Occupant  Type 

Two  columns  for  recording  type  corresponding  to  Uniform  Building  Code 
occupancy  type  list. 


Example:  Student  Dormitory 


H 


Coding  list  consists  of  the  following  types  of  land  use  and  also 
the  standard  building  code  contained  in  this  appendix  section. 

Types  of  Land  Use 


1-F 

one  family 

FIR 

fire  station 

2-F 

two  family 

SEW 

sewage  treatment  plant 

3-F 

three  family 

RGS 

regulation  station 

M-F 

multiple  family 

STD 

students 

SOR 

sorority 

LIB 

library 

FRT 

fraternity 

PRK 

park 

COM 

commercial 

OFF 

office 

CHH 

CHR 

church 

NUR 

nursing  home 

GAR 

garage 

POF 

post  office 

SCH 

school 

WTP 

water  plant 

APT 

apartment 

POL 

state  police  post 

LOT 

parking  lot 

MOT 

motel 

COP 

co-op 

UTL 

utility 

HAL 

hall 

4.  Age  of  Structure 

Two  columns  to  record  year  of  original  construction 


Example:  1948  ] 4[  8 | 


5.  Square  Footage 

Square  feet  of  buildings  will  be  designated  in  hundreds  and  will 
require  four  columns. 


Example:  1200  sq.  ft. 


1 2 


The  number  of  square  feet  given  in  each  case  is  only  for  the  ground  floor. 
To  get  total  square  footage  of  a structure,  square  feet  must  be  multiplied 
by  the  number  of  stories. 
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6.  Number  of  Stories 


Two  columns  to  record  actual  height. 

Example:  17  story  building 


1 7 


Notice  that  the  number  of  stories  is  broken  down  into  fractions , 
i.e.,  1 1/4,  1 1/2,  1 3/4,  etc.  Bi-level  homes  account  for  this 
breakdown.  For  purposes  of  this  study  anything  over  1 story,  such 
as  1 1/2  stories,  would  be  considered  2 stories. 


7.  Structural  Rating 
Requires  one  column. 

Example:  Single  family  home 

8.  Number  of  Persons  Endangered 

Four  columns  recording  actual  number  of  persons  in  structure, 
Example:  3 persons  | | [3 


9 .  Node 


Five  columns  recording  node  will  be  recorded  with  the  first  column 
designating  the  political  jurisdiction  using  the  following  key: 


A = 1 

East  Lansing  City  Limits 

L = 2 

City  of  Lansing 

M = 3 

Meridian  Township 

T = 4 

Lansing  Township 

U = 5 

University  Property  in  East  Lansing 

X = 6 

University  Property  outside  of  East  Lansing 

Example:  M 0003 

3 |0  0 |0  |3 
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GROUP 

A 

B 


C 

D 


DESCRIPTION  OF  OCCUPANCY 


Any  assembly  building  with  a stage  and  an  occupant  load  of 
1000  or  more  in  the  building. 

1- Any  assembly  building  with  a stage  and  an  occupant  load  of 
less  than  1000  in  the  building. 

2 - Any  assembly  building  without  a stage  and  having  an  occupant 
load  of  300  or  more  in  the  building 

3 - Any  assembly  building  without  a stage  and  having  an  occupant 
load  of  less  than  300  in  the  building,  including  such  buildings 
used  for  school  purposes  less  than  eight  hours  per  week. 

4- Stadiums,  reviewing  stands,  and  amusement  park  structures  not 
included  within  Group  A nor  Divisions  1,  2,  and  3,  Group  B 
Occupancies . 

Any  building  used  for  school  or  day-care  purposes  more  than 
eight  hours  per  week,  involving  assemblage  for  instruction, 
education,  or  recreation,  and  not  classified  in  Group  A or 
Divisions  1 and  2,  Group  B,  Occupancies. 

1- Mental  hospitals,  mental  sanitariums,  jails,  prisons,  reforma- 
tories, houses  of  correction,  and  buildings  where  personal 
liberties  of  inmates  are  similarly  restrained. 

2- Nurseries  for  full-time  care  of  children  under  kindergarten 
age.  Hospitals,  sanitariums,  nursing  homes  with  nonambulatory 
patients,  and  similar  buildings  (each  accommodating  more  than 
five  persons) . 

3- Nursing  homes  for  ambulatory  patients,  homes  for  children  of 
kindergarten  age  or  over  (each  accommodating  more  than  five 
persons) . 

1- Storage  and  handling  of  hazardous  and  highly  inflammable  or 
explosive  materials  other  than  flammable  liquids. 

2- Storage  and  handling  of  Class  I,  II,  and  III  flammable  liquids, 
as  specified  in  U.B.C.  Standard  No.  9-1-64;  dry  cleaning  plants 
using  flammable  liquids,  paint  stores  with  bulk  handling;  paint 
shops  and  spray  painting  rooms  and  shops. 
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GROUP 


F 


G 


H 


DESCRIPTION  OF  OCCUPANCY 


3- Woodworking  establishments,  planning  mills  and  box  factories; 
shops,  factories  where  loose,  combustible  fibers  or  dust  are 
manufactured,  processed,  or  generated;  warehouses  where  highly 
combustible  material  is  stored. 

4 - Repair  garages. 

5 - Aircraft  repair  hangar 

1- Gasoline  and  service  stations,  storage  garages  where  no  repair 
work  is  done  except  exchange  of  parts  and  maintenance  requiring 
no  open  flame,  welding,  or  the  use  of  highly  flammable  liquids. 

2- Wholesale  and  retail  stores,  office  buildings,  drinking  and 
dining  establishments  having  an  occupant  load  of  less  than 
100,  printing  plants,  municipal  police  and  fire  station, 
factories  and  workshops  using  material  not  highly  flammable  or 
combustible,  storage  and  sales  rooms  for  combustible  goods, 
paint  stores  without  bulk  handling. 

3- Aircraft  hangars  where  no  repair  work  is  done  except  exchange 
of  parts  and  maintenance  requiring  no  open  flame,  welding,  or 
the  use  of  highly  flammable  liquids. 

Ice  plants,  power  plants,  pumping  plants,  cold  storage  and 
creameries,  factories  and  workshops  using  incombustible  and 
nonexplosive  materials.  Storage  and  sales  rooms  of  incombustible 
and  nonexplosive  materials. 

Hotels  and  apartment  houses.  Convents,  monasteries  (each 
accommodating  more  than  10  persons) . 

1 - Apartments 

2- Frat.,  sor.,  coop. 

3- Student  dormitories 

4- Hotels,  motels 
Dwellings  and  lodging  houses. 

1-Single  family 
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GROUP 


J 


DESCRIPTION  OF  OCCUPANCY 


2 - Rooming  house 

3 - Duplex 

1 - Private  garages,  sheds  and  minor  buildings  used  as  accessories 
only  when  not  over  1000  square  feet  in  area. 

2- Fences  over  6 feet  high,  tanks  and  towers. 


(From  Uniform  Bldg.  Code,  1964  Edition  Pg.  46-48,  Table  No.  5-A.) 


Appendix  E - Data  Tape  Programs 


c 

THIS  PROGRAM  CREATES  NETWORK  TAPE 

A 

05 

c 

THIS  PROGRAM  ENUMERATES  ALL  LINKS 

A 

10 

COMMON  NUDE  A ( 3000  ) * I NODE  (3000)  , NODEtS  ( 3000  ) ,JNGDE(3000>  • D 1ST ( 3 00  0 ) * 

A 

20 

1DIS(3000) , N E X T (30C0) »LOC( 1500) *IDIC( 1500) 

A 

30 

INTEGER  TAP  , R » P 

A 

40 

R = 5 

A 

50 

P = 6 

A 

60 

-T  AP  = 7 

A 

70 

L = 0 

A 

80 

110 

L = L+1 

A 

90 

READ  (R,280,END=115)  INODE(L)  ,JNODE(L)  ,DIS(L) 

A 

95 

K = L + 1 

A 

100 

INODEt  K) =JNODE ( L) 

A 

110 

JNODEt  K)  = INODE  ( L )< 

A 

120 

DIS(K)=DIS(L) 

A 

130 

L =K 

A 

140 

GO  TO  110 

A 

150 

c 

THIS  PROGRAM  SORTS  ALL  LINKS 

A 

160 

115 

I F I RST= 1 

A 

170 

NLINK=L-1 

A 

180 

LAST  = 1 

A 

190 

DO  120  1=1, L 

A 

200 

120 

NEXT ( I ) =0 

A 

210 

DO  200  I=2»NLINK 

A 

220 

I COUNT  = 0 

A 

230 

M=I  FIRST 

A 

240 

130 

IF  ( INODEt  I ) - I NODE ( M ) ) 170,140,180 

A 

250 

1 AO 

IF  (JNODEt I )-JNODE(M) ) 170,150,180 

A 

260 

150 

IF  (M.NE.LAST)  NEXT ( I ) =NEXT ( M ) 

A 

270 

160 

I F ( M.EQ.LAST ) LAST=  I 

A 

280 

NEXT ( M ) = I 

A 

290 

GO  TO  200 

A 

300 

170 

IF  (M.NE. IFIRST)  NEXT(J)=I 

A 

310 

IF  (M.EQ. IFIRST)  IFIRST=I 

A 

320 

NEXT ( I ) =M 

A 

330 

GO  TO  200 

A 

340 

180 

IF  (M.EQ.LAST)  GO  TL  160 

A 

350 

I COUNT  = I CO'UN  T + 1 

A 

360 

IF  ( ICOUNT.LE. I ) GO  TO  190 

A 

370 

WRITE  ( P , 3 5 0 ) I ,M» IFIRST  , LAST  , (NEXT (J)  ,J=1  , I ) 

A 

380 

STOP 

A 

390 

190 

J = M 

A 

400 

M=NEXT (M ) 

A 

410 

GO  TO  130 

A 

420 

200 

CONTINUE 

A 

430 

M=IFIRST 

A 

440 

DO  210  I = 1 , N L I NK 

A 

450 

NODEA ( I ) = I NODE ( M ) 

A 

460 

NODEB ( I ) = J N 0 D E ( M ) 

A 

470 

D I S T ( I ) = D I S ( M ) 

A 

480 

M = N E X T (M) 

A 

490 

210 

CONTINUE 

A 

500 

C 

THIS  PROGRAM  RENUMBERS  ALL  LINKS 

A 

510 

(1 

A 

520 

LONODE=NODEA ( 1 ) 

A 

530 

IDICt 1 ) = LONODE 

A 

540 

LOCt  1 ) =1 

A 

550 

NODEA ( 1 ) =1 

A 

560 
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DO  230  1=2  » N L I N K 

A 

570 

K = K+1 

A 

580 

IF  (NODEA( I ) .EQ.LONODE)  GO  TO  220 

A 

590 

LONODE=NODEA  ( I ) 

A. 

600 

NODEAf I ) =K 

A 

610 

I D I C ( K ) =LONODE 

A 

620 

LOC ( K ) = I 

A 

630 

GO  TO  230 

A 

640 

220 

K = K-  1 

A 

650 

NODEAf I ) =K 

A 

660 

230 

CONTINUE 

A 

670 

DO  260  1 = 1 * NL I NK 

A 

680 

DO  240  J = 1 » K 

A 

690 

IF  (NODEBf I ) .EQ.IDIC(J) ) GO  TO  250 

A 

700 

240 

CONTINUE 

A 

710 

250 

NODEBf  I ) =J 

A 

720 

260 

CONTINUE 

A 

730 

N N 0 D E = K 

A 

740 

WRITE  ( P * 3 6 0 ) 

A 

750 

WRITE  ( P » 370 ) ( IDICf  I ) , 1 = 1, K) 

A 

760 

WRITE  ( P » 320 ) 

A 

770 

i 4 

KK=K+1 

A 

780 

LOC ( KK) =NLI NK+1 

A 

790 

WRITE  ( P * 3 7 0 ) ( LOC ( I ) * I = 1 » KK ) 

A 

800 

WRITE  (TAP)  ( LOC ( I ) * I = 1 * KK ) 

A 

810 

WRITE  (TAP)  ( IDICf  I ) » I = 1 » K ) 

A 

820 

WRITE  ( P » 330 ) 

A 

830 

WRITE  ( P » 340 ) 

A 

840 

DO  270  1 = 1 > N L I NK 

A 

850 

I B = NODEB ( I ) 

A 

860 

I C=NODEA { I ) 

A 

870 

WRITE  ( P , 290  ) IDICf  IC)  , IDICf  IB) ,DIST(  I ) 

A 

880 

WRITE  (TAP)  NODEAf I ) .NODEBf I ) ,DIST( I ) 

A 

890 

270 

CONTINUE 

A 

900 

END  FILE  7 

A 

910 

WRITE  ( P , 300 ) NLINK 

A 

920 

WRITE  ( P » 3 1 0 ) NNODE 

A 

930 

STOP 

A 

940 

C 

A 

950 

280 

FORMAT  ( 29X  , I4,2X , I4,32X, F5. 2 ) 

A 

960 

290 

FORMAT  ( IX,  I 5,2X, I5,6X,F5.2) 

A 

970 

300 

FORMAT  (1H022HTHE  NUMBER  OF  LINKS  I S I 5 ) 

A 

980 

310 

FORMAT  ( 1H022HTHE  NUMBER  OF  NODES  I S I 5 ) 

A 

990 

320 

FORMAT  ( 1H1 , 23HVALUES  OF  LOC  TO  FOLLOW) 

A 1000 

330 

FORMAT  ( 1H1,41HLISTING  OF  RENUMBERED  LINK  DATA 

TO  FOLLOW) 

A1010 

340 

FORMAT  ( 1H0 ,2X,6HNODE  A,2X,6HNODE  B,2X,10HLINK 

DIST. ) 

A 1 0 2 0 

350 

FORMAT  ( 22H  ERROR  IN  SORTING-I  = , I 4 , 9H I F I RS T = 

* , 14, 8H. LAST  = ,14 

A 1 0 3 0 

1/35H  THE  VALUES  OF  NEXT  ARE  AS  FOLLOWS-/ ( 2 5 I 5 ) ) 

A 1040 

360 

FORMAT  ( 1H1  , 30HVALUES  OF  DICTIONARY  TO  FOLLOW) 

A 1 0 5 0 

370 

FORMAT  (1H0,20I6) 

A 1060 

END 

A1070- 
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C THIS  PROGRAM  CREATES  STRUCTURE  INDEX  B 05 

DIMENSION  DUM ( 9 ) B 10 

INTEGER  P,R,TAP  B 20 

R = 5 B 30 

P=6  B AO 

T AP  = 7 B 50 

I COUNT  = 0 B 60 

110  READ  ( R , 1 40  » END= 1 1 5 ) (DUM( 1)  *DUM(2) , JDUM » (DUM(J) » J = 4 » 8 ) ) » I OLD  B 65 

WRITE  (TAP.1A0)  ( DUM ( 1 ) , DUM ( 2 ) * JDUM * ( DUM ( J ) , J = A , 8 ) ) , I OLD  B 70 

I COUNT= ICOUNT+1  B 80 

GO  TO  110  B 90 

115  END  FILE  7 B 100 

WRITE  ( P » 1 3 0 ) ICOUNT  B 110 

WRITE  (P.120)  B 120 

STOP  B 130 

C • B 1 AO 

120  FORMAT  (20X.26HTHE  JOB  HAS  BEEN  DONE  WELL)  B 150 

130  FORMAT  (25H  THE  VALVE  OF  ICOUNT  IS  ,15)  B 160 

1 AO  FORMAT  ( F7 . 1 , 1 X , F5 . 1 , AX  , AA  , 3X  , F A . 0 , 2X  , F5  . 1 , 3X  , F A . 1 , AX  » F 3 . 0 IX  , F6 . 0 B 170 
1.19X.I5)  B 180 

END  B 190- 
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Appendix  F - Fire  Station  Location  Program 


c 


1 10 


120 


130 

140 

150 

C 

160 


170 

180 

190 

200 

210 


THIS  PROGRAM  EVALUATES  SITES  FOR  THE  LOCATION  OF  A FIRE  STATION.  A 10 
COMMON  MAX , NNODE ,NL I NK , I RROR  , IOPT 1 , I OPT2 »NUNI TS ♦ JDUM, I TAPE , IOLD .NO  A 20 
1 D E ( 1500)  * TM  1 ( 1500 ) , DIST ( 1500  ) *XX ( 500 ,9  ) ,DUM ( 9 ) *IC0DE(70)  * LOC ( 500 ) » A 30 
2IDICT (500)  > X ( 500  » 9 ) » A ( 9 ) * F ( 5 0 0 ) ,TOTF(500)  * DETROY ( 10),TRITI(500)*KP  A 40 
3ATH ( 5 00  ) * SEQNOD ( 150)  *TOTG(500)  *LIMIT( 10*10)  » NPA  TH ( 150  > » NBACK ( 150  ) ♦ A 5 0 
4ICK(10),ICHO(10),TITLE(12),IPOL(500),DISTA(50Q),Y<500)»YY(500>»ASS  A 60 
5IGN ( 70  ) , TOTDIS*SUMT IM,  IDEST(10),NSTR(3), ISUM(  10.10)»JSID(500)i  I MAX  A 70 


INTEGER  P*R, FUTURE,  DETROY,  TAP  A 80 

R = 5 A 85 

P = 6 A 86 

READ  ( R » 4 1 0 ) NRUNS, IOPT1 , IOPT2.FUTURE  A 90 

READ  (R»420)  MAX , NNODE , NL I NK  , TAP  A 100 

READ  ( R » 47  0 ) ( I CODE ( I ) » I = 1 » 64 ) A 110 

READ  ( R > 480 ) (ASSIGN!  I ) * 1 = 1,64)  A 120 

READ  ( R » 490 ) (A(I),I=1,9)  A'l30 

WRITE  ( P » 5 0 0 ) A 140 

DO  110  1=1,60  A 150 

WRITE  ( P 9 5 1 0 ) I CODE ( I ) , ASS  I GN ( I ) A 160 

CONTINUE  A 170 

DO  120  1=1, NNODE  A 180 

Y ( I ) =0 .0  A 190 

DO  120  J = 1 , 9 A 200 

X ( I , J ) =0.0  A 210 

ND  = NNODE+ 1 A 220 

READ  (TAP)  (LOC(KO)  ,KO=l»ND)  A 230 

RFAD  (TAP)  ( IDICT (KK)  ,KK  = 1 , NNODE)  A 240 

I 7 A P E = 0 A 250 

NUN  I TS= 1 A 260 

IF  ( FUTURE. NE. 1 ) GO  TO  150  A 270 

DO  130  1=1,10  A 280 

DO  130  J = 1 » 3 A 290 

I S U M ( I , J ) =0  A 300 

READ  ( R ,450  ) (NSTR(  I ) , 1 = 1,3)  A 310 

DO  140  1=1,10  A 320 

READ  ( R , 460 ) ( DETROY ( I ) , ( L I M I T ( I , J ) , J = 1 , 3 ) ) A 330 

NR=0  A 340 

NR=NR+ 1 A 350 

READ  ( TAP,660,END=160)  NODA , NODE ( NR ) , D I ST ( NR ) , TM 1 ( NR ) A 355 

GO  TO  150  A 360 

THIS  SET  OF  INSTRUCTIONS  DELETES  STRUCTURES  FROM  CONSIDERATION  A 370 

READ  ( TAP,640,END  = 235)  ( DUM  ( 1 ) , DUM  ( 2 ) , JDUM , ( DUM ( I DUM ) , I DUM  = 3 , 7 ) ) * A 375 

HOLD  A 376 

IF  ( FUTURE. NE. 1 ) GO  TO  210  A 380 

DO  170  I D ~ 1 , 10  A 390 

IF  (DETROY! ID) .EQ. IOLD)  GO  TO  180  A 400 

CONTINUE  A 410 

GO  TO  210  A 420 

DO  190  J J= 1 ,3  A 430 

IF  ( JDUM . EQ . NSTR ( J J ) ) GO  TO  200  A 440 

CONTINUE  A 450 

GO  TO  210  A 460 

I S U M ( ID»JJ)  = ISUM( I D , J J ) + 1 A 470 

IF  (LIMIT ( I D » JJ ) .LT. ISUM( I D » JJ ) ) GO  TO  210  A 480 


GO  TO  160  A 490 
DUM ( 8 ) = 1 . 0 A 500 
IF  ( DUM ( 3 ) .LT.80  ) GO  TO  220  A 510 
DUM  ( 3 ) =180-DUM  ( 3 ) 'A  520 
GO  TO  230  F-1  A 530 


220 

230 


235 
24  0 


250 


260 


270 


280 


290 


300 

310 

320 


330 

340 

350 

360 


DUM ( 3 ) =80-DUM ( 3 ) A 540 

I TAPE= I TAPE  + 1 A 550 

CALL  HOUSE  A 560 

GO  TO  160  A 570 

IF  ( FUTURE. NE.l ) GO  TO  250  A 580 

READ  ( R , 700  » END=2  5 0 ) ( DUM ( 1 ) » NUN  I TS » JDUM * ( DUM ( I DUM ) » I DUM= 3 » 7 ) ) » A 585 

HOLD  A 586 

DUM ( 8 ) =1 .0  A 590 

I TAPE= ITAPE+1  A 600 

CALL  HOUSE  A 610 

GO  TO  240  A 620 

WRITE  (P.580)  A 630 

WRITE  ( P 5 540 ) NNODE, NLINK  A 640 

..RITE  ( P * 600  ) ITAPE  A 650 

WRITE  (P.650)  A 660 

WRITE  ( P , 6 70 ) A 670 

DO  260  I=l,NNODE  A 680 

WRITE  (P.680)  I . ( X ( I . J ) . J=1 ,9 ) A- 690 

CONTINUE  A 700 

WRITE  (P.690)  (A(I). 1=1.9)  A 710 

DO  360  NSTATN=1 .NRUNS  A 720 

DO  270  1 = 1 , NNODE  A 730 

TO TG ( I ) =0 . 0 A 740 

TOTF ( I ) =0.0  A 750 

F ( I ) =0.0  A 760 

DO  280  NK= 1 » NNODE  A 770 

YY(NK)=Y(NK)  A 780 

DO  280  K N = 1 ,9  A 790 

XX (NK.KN) =X( NK.KN)  A 800 

READ  (R.440)  ( I MAX , T I T L E ( I ) , I = 1 , 1 2 ) A 810 

WRITE  (P.430)  ( TITLE ( I ) , 1=1 , 12 ) A 820 

WRITE  (P.580)  A 830 

READ  (R.630)  ( I CK ( I B ) , I B= 1 , I MAX ) A 840 

READ  (R.610)  ( I CHO ( IK) , IK=1,IMAX)  A 850 

READ  < R , 630 ) ( I DEST ( I D ) » I D= 1 , I MAX ) A 860 

WRITE  (P.570)  A 870 

DO  290  I R I G I N= 1 » I MAX  A 880 

WRITE  (P.590)  I R I G I N » I CK ( IRIGIN)  A 890 

READ  (R.560)  JSS  A 900 

IF  (JSS.EQ.O)  GO  TO  320  A 910 

READ  (R.620)  ( JS I D ( JS ) , JS=  1 , NNODE ) A 920 

D(  310  JS=1,  NNODE  A 930 

IF  (JSID(JS) .EQ.l)  GO  TO  310  A 940 

YY ( JS ) =0.0  A 950 

DO  300  1=1,9  A 960 

XX( JS.I  ) =0 . 0 A 970 

CONTINUE  * A 980 

READ  (R.560)  NPOL  A 990 

IF  (NPOL.EQ.O)  GO  TO  350  A1000 

READ  ( R , 550 ) POLICY  A1010 

DO  340  INEW=1, NNODE  A1020 

READ  (R.620)  IPOL(INEW)  A1030 

IF  ( IPOL(  INEW)  .EQ.O)  GO  TO  340  A1040 

DO  330  I DUM= 1,9  A1050 


XX ( I NEW  , I DUM )=POLICY*XX(  INEW,  I DUM ) 

CONTINUE 

CALL  CARTIM 


A 1060 
A 1 070 
A 1080 


CONTINUE 
WRITE  (P.580) 

I F ( IOPT2.EQ. 1 ) GO  TO  370 
GO  TO  400 


A1090 
A1 100 
A 1 1 1 0 
A1120 
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370 

WRITE  i 

[ P ,520  ) 

A 1 1 3 0 

DO  380 

I DT= 1 , NNODE 

A1 140 

DO  380 

I X=  1 ,8 

A 1 1 5 0 

F<  IDT  ) = 

= F ( IDT  ) + (X(  IDT, IX  ) * A ( IX)  ) 

A1 160 

380 

CONTINUE 

A1 170 

DO  390 

1=1  , NNODE. 

A 1 1 8 0 

390 

WRITE  ( 

:P»530  ) I , I D I C T ( I ) ,F(  I ) 

A 1 1 90 

400 

STOP 

A 1 200 

C 

A 1 2 1 0 

410 

FORMAT 

( 13,3 ( 2 X , I 1 ) ) 

A 1 2 2 0 

420 

FORMAT 

(315,12) 

A 1 2 30 

430 

FORMAT 

( 1H1 » 30X , 12A6  ) 

A 1 240 

440 

FORMAT 

( I 5 , 1 OX , 12A6 ) 

A 1 2 5 0 

450 

FORMAT 

( 3A4  ) 

A 1 2 60 

460, 

FORMAT 

( 16, IX, 313) 

A 1 270 

470 

FORMAT 

( 1 6 ( A4 ) ) 

A1280 

480 

FORMAT 

( 16F4.2 ) 

A 1 290 

490 

FORMAT 

(9F6.3) 

A 1 3 00 

500 

FORMAT 

( 20H  LAND  USE  VALUE) 

A 1 3 1 0 

510 

FORMAT 

( 5X,A4,5X,F4.2 ) 

A 1 320 

5 20 

FORMAT 

( 2X  » 1 OHD I CT IONARY , 5X, 16HWEIGHT ING  FACTOR ) 

A 1 3 30 

530  - 

FORMAT 

<2X, I3»3X,I5»5X,F10.2) 

A 1 340 

540 

FORMAT 

( 10X , 19HTHE  NO.  OF  NODES  IS,I6»25H,  AND  THE  NO.  OF  LINKS  IS 

A 1 3 50 

1.17) 

A 1 3 60 

550 

FORMAT 

(F4.2  ) 

A 1 3 70 

560 

FORMAT 

( 13  ) 

A 1 3 80 

570 

FORMAT 

( 36H  NODES  CHOSEN  FOR  POSSIBLE  LOCATIONS) 

A 1 3 90 

580 

FORMAT 

( 1 HO , ///  ) 

A 1 400 

590 

FORMAT 

( 1H  , I 5 , IX, 1H ( , I 5 , 1H)  ) 

A 1 4 1 0 

600 

FORMAT 

( 1 1 H THERE  ARE ,16, 28H  STRUCTURES  BEING  CONSIDERED) 

A 1 420 

610 

- FORMAT 

( 1 0 ( 4 X , I 1 ) ) 

A 1430 

620 

FORMAT 

(8011) 

A 1 440 

630 

FORMAT 

( 10 ( IX, 14) ) 

A 1 45  0 

640 

FORMAT 

(F7.1,1X,F5.1»4X,A4,3X,F4.0,2X,F5.1,3X,F4.1,4X,F3.0»1X»F6.0 

A 1 460 

1 » 19X  ,15) 

A 1 47  0 

650 

FORMAT 

( 1H1 ,40X,16HPARAMETER  VALUES) 

A 1 48  0 

660 

FORMAT 

( 2 I 5 , 2F4. 2 ) 

A1490 

670 

FORMAT 

( 1 IX . 44HSTR  VAL  LAND  VAL  STR  AGE  SQ  FT  STORIES, 40H 

A1500 

1RATING 

PEOPLE  CONSTANT  LAND  USE) 

A 1 5 1 0 

680 

FORMAT 

(2X,I4,2X,8(F9.0)»5X,F6.2) 

A 1 52  0 

690 

FORMAT 

( 8H  WEIGHTS ,3X,9 ( F6.3 ,3X  ) ) 

A 1 5 3 0 

700 

FORMAT 

(F7.1,2X,I6»4X,A4,2X,F5.1,2X»F5.1»3X,F4.1,4X,F3.1,1X,F4.0»1 

A 1 540 

19X.I5  ) 

. A 1 550 

END 

A 1 560 

SUBROUTINE  CARTIM 

C 10 

C 

THIS  SUBROUTINE  IS  THE  SHORTEST  PATH  ALGORITHM 

C 20 

COMMON 

MAX, NNODE.NL INK, IRROR , I OP T 1 , I OP T 2 , NUN  I T S , J DUM , ITAPE.IOLD.NO 

C 30 

1DE( 1500  ) »TM1 ( 1500 ) »DIST ( 1500 ) ,XX ( 500,9 ) »DUM( 9)  . I CODE ( 70 ) , LOC ( 5 00 ) , 

C 40 

2IDICT ( 500 ) ,X ( 500,9 ) ,A ( 9 ) ,F( 500 ) ,TOTF( 500 ) .DETROY ( 10 ) ,TRITI { 500 ) ,KP 

C 50 

3ATH ( 500  ) , SEQNOD ( 150)  *TOTG(500) .LIMIT ( 10,10)  »NPATH(150*  » NBACK ( 150) , 

C 60 

41 CK ( 10  ) 

, I CHO ( 10)  , TITLE!  12)  ,IPOL(500)  ,DISTA(500)  * Y ( 5 00  ) »YY(500)  »ASS 

C 70 

5IGN( 70  ) 

»TOTDIS»SUMTIM,  IDEST(10),NSTR(3), ISUM(  10,l0)»JSID(500)  , I MAX 

C 80 

INTEGER 

SEQNOD,  ANODE, P,R, TAP 

C 90 

N«NNODE 

C 100 

DO  130 

I A=1 , IMAX 

c no 

DO  110 

I C= 1 , NNODE 

C 120 

IF  ( I C K ( I A ) . EQ -IDICT(IC)  ) GO  TO  120 

C 130 

110 

CONTINUE 

C 140 

WRITE  ( P , 560 ) IA 

C 150 

GO  TO  130 

C 160 

120 

I CK ( I A) 

= IDICT ( IC ) 

C 170 
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130  CONTINUE  C 180 

DO  410  10=1 , IMAX  C 190 

IRIGIN=ICK( 10)  C 200 

DO  140  J=1 .NNODE  C 210 

TRI T I ( J ) =9000.  C 220 

140  KP ATH ( J ) =0  C 230 

TRITI ( IRIGIN)=0.  C 240 

DO  130  MM= 1 » MAX  C 250 

SEQNOD ( MM ) = 0 C 260 

150  NPATH ( MM ) =0  C 270 

NODEA=IRIGIN  C 280 

I PUT= 1 C 290 

I GET=  1 C 300 

I COUNT  = 0 C 310 

AT  I ME  = 0 • 0 C 320 

C THIS  SET  OF  INSTRUCTIONS  SEARCHES  ALL  LINKS  LEADING  FROM  A NODE.  C 330 

160  NBR=LOC ( NODEA ) C 340 

NEXT=LOC ( NODEA+1 ) -1  C 350 

DO  180  J J = NBR » NEXT  C 360 

IF  (NODE(JJ) .EQ. IRIGIN)  GO  TO  180  C 370 

NODEB=NODE ( J J ) C 380 

TIME=ATIME+TM1< JJ)  C 390 

IF  (TIME. GE. TRITI (NODEB) ) GO  TO  180  C 400 

TRITI (NODEB)=TIME  C 410 

KPATH(NODEB) =NODEA  C 420 

SEQNOD ( IPUT ) =NODEB  C 430 

I COUNT= I COUNT+ 1 C 440 

IF  ( ICOUNT.LE.MAX ) GO  TO  170  C 450 

WRITE  (P.520)  ICOUNT  C 460 

GO  TO  410  C 470 

170  IPUT=IPUT+1  C 480 

IF  ( IPUT.EQ.MAX+1 ) IPUT=1  C 490 

180  CONTINUE  C 500 

IF  ( ICOUNT. EQ.O ) GO  TO  190  C 510 

NODEA  = SEQNOD( I GET)  C 520 

ATIME=TRITI (NODEA)  C 530 

ICOUNT=ICOUNT-l  C 540 

IGET=IGET+1  C 550 

IF  ( IGET.EQ.MAX+1 ) IGET=1  C 560 

GO  TO  160  C 570 

190  SUMT IM=0.0  C 580 

DO  200  J = 1 . NNODE  C 590 

IF  ( TRITI  (J)  .GE. 9000.  ) GO  TO  200  C 600 

IF  (J.EQ. IRIGIN)  GO  TO  200  C 610 

SUMTIM=SUMTIM+ TRITI (J)  C 620 

200  CONTINUE  C 630 

DO  320  J=l» NNODE  C 640 

IF  ( IRIGIN. EQ.J)  GO  TO  320  C 650 

I B = MAX  C 660 

MPATH= J C 670 

NP ATH ( IB ) =J  C 680 

210  I B = I B — 1 C 690 

IF ' ( IB.GT.O)  GO  TO  220  C 700 

I RROR=274  C 710 

WRITE  (P.500)  IRROR  C 720 

GO  TO  320  C 730 

220  NPATHf IB) =KPATH ( MPATH ) C 740 

IF  (NPATH( IB). EQ. IRIGIN)  GO  TO  230  C 750 

MPATH=NPATH ( IB ) C 760 

GO  TO  210  C 770 

230  IF  ( ICHO(  10)  .EQ.O)  GO  TO  280  C 780 
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DO  240  I DD= 1 » IMAX  C 790 

IF  ( IDICT ( J ) .EQ. IDEST ( I DD ) ) GO  TO  250  C 800 

240  CONTINUE  C 810 

GO  TO  280  C 820 

250  DO  260  LQ=1,MAX  C 830 

260  NBACK(LQ)=0  C 840 

DO  270  LK= I B .MAX  C 850 

JK=NPATH ( LK ) C 860 

INT= IDICTt JK  ) C 870 

270  NBACK(LK)=INT  C 880 

WRITE  (P.580)  IRIGIN.J  C 890 

WRITE  (P.570)  (NBACK(JQ) »JQ=IB»MAX)  C 900 

280  IF  (IOPT1.EQ.O)  GO  TO  380  C 910 

D I S T A ( J ) =0 . C 920 

T0TDIS=0.0  C 930 

NMAX=MAX~1  C 940 

DO  310  I X= I B . NMAX  C 950 

ANODE=NPATH ( IX ) C 960 

NTRY=LOC ( ANODE ) C 970 

LIMIT1=L0C( ANODE+1 )-l  C 980 

DO  290  JQ  = N TRY  » L I M I T 1 C 990 

IF  ( NPATH ( IX+1 ) . EQ .NODE ( JQ ) ) GO  TO  300  C1000 

290  CONTINUE  C1010 

I RROR=277  C1020 

WRITE  (P.,500)  IRROR  C1030 

GO  TO  310  C 1 040 

300  DISTA(JQ)=DIST(JQ)+DISTA(J)  C1050 

310  CONTINUE  C1060 

320  CONTINUE  C1070 

DO  330  J = 1 . NNODE  C1080 

330  TOTDIS=TOTDIS+DISTA( J)  C1090 

IRGN=IDICT( IRIGIN)  C1100 

WRITE  IP. 480)  C1110 

WRITE  (P.440)  IRIGIN. IRGN  0120 

JR  = 1 C 1 1 30 

JM=  1 5 C 1 1 40 

340  WRITE  (P.450)  ( T R I T I ( JB  ) , JB  = JR  , JM  ) 0150 

JP  = JR+1  5 0160 

JOJM+15  0170 

IF  (N-JM)  350.350.340  0180 

350  WRITE  (P.450)  (TRITI  (JB)  >JB  = JR»NNODE)  0190 

WRITE  (P.540)  IRIGIN. IDICTt  IRIGIN)  .SUMTIM  0 200 

WRITE  (P.470)  0210 

WRITE  (P.460)  IRIGIN, IRGN  0220 

JR  = 1 0230 

JM  = 1 5 0240 

360  WRITE  (P.450)  ( D I S T A ( MB ) , MB  = JR  , JM ) 0250 

JR  = JR+15  0260 

JM  = JM+ 1 5 0270 

IF  (N-JM)  370,370.360  0280 

370  WRITE  (P.450)  ( D I S T A ( MB  ) , MB  = J R , NNODE  ) 0 290 

WRITE  (P.550)  IRIGIN,  IDICTt  IRIGIN)  .TOTDIS  0 300 

380  TOTF(  I R I G I N ) =0  • 0 0310 

TOTG(  I R I G I N ) =0  • 0 0320 

DO  390  IDT=1, NNODE  0330 

F ( I D T ) = 0 . 0 0340 

DO  390  I X=  1 , 8 0350 

390  F ( IDT  )=F  ( IDl  ) + ( X(  IDT  , IX  )*A(  IX)  ) 0360 

DO  400  IDT=1, NNODE  0 370 

IF  (TRITI  ( IDT)  .GE. 9000.0)  GO  TO  400  0380 

TOTF(IRIGIN)=TOTF(IRIGIN)  + (F(IDT)*TRITI(IDT))  0 390 
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400 

410 

420 

430 

C 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 

540 

550 

560 

570 

580 

110 

120 

130 


TOTG(  IRIGI.  N)=TOTG(  IRIGIN)  + (YY(  IDT  ) * T R I T I ( IDT)  ) C1400 

CONTINUE  C1410 

CONTINUE  C 1420 

WRITE  ( P * 5 1 0 ) C 1430 

WRITE  ( P » 490 ) C1440 

DO  430  I = 1 * NNODE  ' C1450 

IF  (ICK(I).EQ.O)  GO  TO  430  C1460 

NEXT= I C 1 470 

NNEXT  = NEXT  + 1 C1480 

DO  420  J=NNEXT *NNODE  C1490 

IF  ( ICK(J) .EQ.O)  GO  TO  420  C1500 

II'  (TOTF(J)  .LT.TOTF(NEXT)  ) NEXT  = J C1510 

CONTINUE  C 1 5 20 

WRITE  ( P » 530  ) NEXT  , IDICT ( NEXT  ) ,TOTF( NEXT ). TOTGf NEXT ) C1530 

TOTF(NEXT)=2.0**35-1  C1540 

CONTINUE  C 1 550 

RETURN  C 1 560 

C 1 570 

FORMAT  ( 40X , 1H0 »26HVALUES  OF  TIME  FROM  OR  I G I N , 1 X , I 5 » 1 X , 1 H ( , I 5 , 1H ) , C1580 

1/)  C 1 590 

FORMAT  (1H  , 6X * 15 ( F5 .2 ,2X  ) ) C1600 

FORMAT  ( 40X , 1H0 , 30HVALUES  OF  DISTANCE  FROM  OR  I G I N , IX , I 5 , IX » 1H  < , I 5 , C1610 

1 1H ) » / ) C1620 

FORMAT  ( 1 HO  » / ) C1630 

FORMAT  ( 1 HO  » //  / ) C1640 

FORMAT  ( 16X > 59HRENUMBERED  ORIGINAL  WEIGHTED  TIMES  SUMMED  C1650 

1 TIMES)  C 1 660 

FORMAT  (42H  ERROR  CONDITION  FOUND  AT  STATEMENT  NUMBER , I 5 , 2 9H . PRO  C1670 
1CESSING  CANNOT  PROCEED.)  C1680 

FORMAT  ( 45H  LISTING  OF  WEIGHTED  SUMS  TO  ALL  DESTINATIONS)  C1690 

FORMAT  ( 46H  THE  SEQUENCE  TABLE  CAPACITY  HAS  BEEN  EXCEEDED/10H  ICOU  C1700 

1 NT  = ,14)  C1710 

FORMAT  (1H  ,15X,I5,8X,I5,3X,F14.2,8X,F14.2)  C1720 

FORMAT  ( 1H0 ,22HTOTAL  TIME  FROM  OR  I G I N , IX , I 5 , 1H ( , I 5 » 1H  ) » 1 X , 2H  I S ■ F 1 0 C1730 

1.2  » 1X»7HMINUTES ) „ C1740 

FORMAT  ( 1H0.26HTOTAL  DISTANCE  FROM  OR  I G I N , 1 X , I 5 , 1 H ( , I 5 » 1 H ) , IX , 2H I S C1750 
1,F10.2»1X,5HMILES)  C1760 

FORMAT  ( 20X , 19HTHE  LOCATION  OF  THE,I3*25HTH  ORIGIN  CANNOT  BE  FOUND  C1770 
1)  C 1 780 

FORMAT  (1H  , 20  ( I 5 , IX ) /1H  , 20 ( I 5 , 1 X ) / 1 H ,20(15, IX))  C1790 

FORMAT  (17H0PATH  FROM  OR  I G I N , I 3 , 1 5H  TO  DEST I NAT  I ON , I 4 ) C1800 

END,  C 1 8 10 

SUBROUTINE  HOUSE  D 10 


COMMON  MAX, NNODE, NLINK, I RROR, IOPT 1 , IOPT2 >NUNI TS, JDUM, I TAPE » IOLD.NO  D 20 
IDE  ( 1500  ) ,TM1  ( 1500  ) ,DIST  ( 15.00  ) ,XX(  500 ,9  ) ,DUM(  9 ) , I CODE  ( 7 0 ) , LOC  ( 5 00  ) , D 3 0 
2IDICT (500)  , X ( 500 , 9 ) »A(9)  , F ( 500 ) ,TOTF(500)  , DETROY ( 10) ,TRITI (500)  .KP  D 40 
3ATH( 5 00  ) ♦ SEQNOD ( 150),TOTG(500)»LIMIT(10,10)  , NPA  TH ( 150)*NBACK(150)»  D 50 
4 1 CK.  ( 10)  , I CHO  ( 10)  , T I TLE  ( 12)  >IPOL(50O)  »DISTA(500)  » Y ( 500  ) *YY(500)  »ASS  D 60 
5 I GN ( 70 ) ,TOTDIS,SUMTIM, IDEST(10),NSTR(3) , ISUM(  10,10),JSID(500),IMAX  D 70 


INTEGER  P D 80 
DO  120‘  1 = 1, NNODE  D 90 
IF  (IDICT(I).EQ.IOLD)  GO  TO  110  D 100 
GO  TO  120  D 110 
I NEW= I D 120 
GO  TO  130  D 130 
CONTINUE  D 140 
WRITE  (6,250)  IOLD  D 150 
WRITE  ( P , 240  ) ITAPE  D 160 
GO  TO  230  D 170 
DO  150  1=1,60  D 180 
IF  ( ICODEt I ) .EQ. JDUM)  GO  TO  140  D 190 
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GO  TO  150 

D 

200 

140 

DUM ( 9 ) =ASS I GN ( I ) 

D 

210 

GO  TO  160 

D 

220 

150 

CONTINUE 

D 

230 

WRITE  ( P * 2 6 0 ) JDUM 

D 

240 

WRITE  ( P » 240 ) ITAPE 

D 

250 

GO  TO  230 

D 

260 

160 

IF  ( DUM ( 5 ) . L T • 3 ) GO  TO  170 

D 

270 

IF  ( DUM ( 5 ) . L T . 5 ) GO  TO  180 

D 

280 

I F ( DUM ( 5 ) .LT.7 ) GO  TO  190 

D 

290 

DUM ( 5 ) =10.0 

D 

300 

GO  TO  200 

D 

310 

170 

DUM ( 5 ) =0.0 

D 

320 

GO  TO  200 

D 

330 

— 3r8  0 

DUM ( 5 ) =2 . 0 

D 

340 

GO  TO  200 

D 

350 

190 

DUM ( 5 ) =4.0 

D 

360 

200 

DO  210  1=1,8 

D 

370 

210 

DUM ( I ) =NUNI TS*DUM( I )*DUM(9) 

' D 

380 

Y ( I N EW ) = 1 . 0 

D 

390 

DO  220  1=1,9 

D 

400 

X( I NEW,  I )=X(  I NEW,  I ) +DUM ( I ) 

D 

410 

220 

CONTINUE 

D 

420 

230 

RL'TURN 

D 

430 

C 

D 

440 

240 

FORMAT  (10X,40HTHE  MISTAKE  OCCURS  ON 

STRUCTURE 

CARD  NO. , 16  ) 

D 

450 

250 

FORMAT  (24H  THE  LOCATION  OF  IOLD,I6 

, 1 7H  CAN 

NOT  BE  FOUND) 

D 

460 

260 

FORMAT  ,( 2 7H  THE  VALUE  FOR  LAND  USE 

' , A4 , 1 8H 

CAN  NOT  BE  FOUND) 

D 

470 

END 

D 

480- 
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